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The zero-range potential is customarily employed in vari- 
ous mean-field calculations of many-body systems in atomic 
and nuclear physics within, correspondingly, Gross-Pitaevskii 
and Skyrme-Hartree-Fock approach. We argue, however, that 
a many-body system with zero-range potentials is unstable 
against clusterization into collapsed three-body subsystems. 
We show that neither the density dependence of the potential 
nor an additional repulsive three-body potential can prevent 
this unexpected correlational collapse if the potentials are of 
zero range. Therefore the zero-range potential can only be 
used in many-body calculations where all three-body correla- 
tions are explicitly excluded. 
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I. INTRODUCTION 

A. Many-body systems with zero-range interaction 

The zero-range potential is a very useful concept for 
physical problems with a distinct separation of scales 
|J. If, say, the typical wavelength of a particle in a 
potential well is much larger than the range of the po- 
tential the latter can be approximated by a zero-range 
(pseudo)potential Vq which acts on the particle's wave- 
function V>(r) as |]j§ 

VW(r)=t 5(r)£r^(r), (1) 
or 

where the strength of the potential to is expressed in 
terms of the mass m and the s-wave scattering length a 
of the particle[] as <o = —4:Trh 2 a/m. For a well-behaving 
ip(t) the operator J^r can be replaced with unity. 

The zero-range potentials are widely used in different 
fields of physics as a practically convenient form of the ef- 
fective interaction. In particular it is routinely employed 
in various Hartree-Fock calculations in atomic |^-|5| , nu- 
clear [p|-[ro||, and astro-physics ]TlJ . The Hamiltonian of 
a many-body system is approximated by an operator^] Q 

N / h 2 \ 1 
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1 kcot(6 ) = l/a + 0(k 2 ). 

2 The full Skyrme potential includes also momentum- 
dependent terms which, however, are not important for the 
collapse discussed here. 



where is the coordinate of the i-th particle and N is 
the total number of particles. If the system is exposed 
to an external field a corresponding term is added to the 
Hamiltonian. 

The Hartree-Fock approximation to the wavefunction 
of a system of bosons leads then to the Gross-Pitaevskii 
energy functional Jl2 ], 

EB= J dr (^ Ts(r) + f °^ (r) ) ' (3) 

where Tg (r) is the kinetic energy density and n g (r) is the 
particle density of the system of bosons. For a system of 
fermions one obtains instead the Skyrme-Hartree-Fock 
energy functional |l5| ] which, taking symmetric nuclear 
matter as an example, can be written in terms of the 
corresponding fermionic densities Tf and np as 

Ef = I dr (£ TF(r) + ^ (r) )- (4) 

The ground state wavefunction of the Hamiltonian 
minimizes (within the Hartree-Fock approximation) the 
corresponding energy functional. 

B. Mean- field collapse of a many-body system with 
attractive zero-range interaction 

The kinetic energy Tp of symmetric nuclear matter can 
be expressed in terms of the Fermi momentum, kp = 
{3ir 2 n F /2) 1/3 , as H| 

2 3 5 3/3vr 2 \ 2/3 5/3 

Inserting this in Eq.(^) shows that for large nj? the po- 
tential energy term, proportional to np, prevails over the 

5/3 

kinetic energy term, proportional to rip . Therefore in 
the case of attraction, to < 0, the degeneracy pressure 
of the fermionic gas can not withstand the attraction 
and the system collapses: the minimum of the energy, 
Eo = — oo, is reached at infinitely high density np — oo. 

Since there is no degeneracy pressure in a system of 
bosons the collapse in this case is even more severe [|l5| . 

The mean field collapse, however, can be easily re- 
moved within the same approach by extending the Hamil- 
tonian (^) with either a repulsive three-body zero-range 
potential 

W 3 = J2 M(r ( - rj-)<5(ri - r fc ) , (6) 

i>j>k 
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or a density dependent zero-range potential 



W 3 



i>j 



i 3 n Q (r i )5(r l - r.,-) , 



(7) 



which results in an additional term in the energy func- 
tional proportional to n 3 or n a+2 . The energy function- 
al © and (|) acquire then a well-defined minimum at a 
finite saturation density (l| . 

This mean-field collapse can also occur for some fi- 
nite range interactions, but an appropriate density de- 
pendence can in this case also remove the collapse flq|. 



/ 1 rrijrnk 
to rrij + mfc 



Oj - r fc ) , 



(8) 



_ _ J 1 mj(mj + m k ) ^ _ rrijTj + m fc r fc ^ 



to rrii + TOj + TOfe ^ " rrij + J 
psm(a t ) = Xi , pcos(a i ) = y t , 

where {i,j, k} is a cyclic permutation of {1,2,3} and to 
is an arbitrary mass scale. The set of angles fli consists 
of the hyper- angle on and the four angles Xj/|xj| and 
y i / 1 y -i I ■ T ne kinetic energy operator T is defined as 



C. Correlational collapse of a many body system 
with zero-range potentials 

Curiously enough the rigorous solution of a three-body 
problem with zero-range potentials also exhibits a col- 
lapse known as the Thomas effect [ JL7| |. However, this is 
a different type of collapse, where infinitely many bound 
states appear with exceedingly large binding energies and 
exceedingly small spatial extension. The many body 
system with Hamiltonian (Q) is therefore also subjected 
to this non-mean-field correlational collapse - clusteriza- 
tion into collapsed three-body subsystems. Rather sur- 
prisingly this undesirable feature of the Hamiltonian ^) 
seems to have been unnoticed so far. 

The additional potential W3 removes the mean-field 
collapse of the many-body system. But unless it also re- 
moves the correlational collapse, the Hamiltonian (g) will 
still have this unpleasant property of not having a finite 
ground state energy. The many on-going investigations 
would then be disturbingly close to a divergence, where 
even a small admixture of three-body correlations might 
influence the solution by a substantial and uncontrollable 
amount. 

We shall show in this letter that neither the three- 
body zero-range potential (^) nor the density dependent 
zero-range potential (0) is actually able to remove the 
collapse of the many-body system into collapsed three- 
body clusters. 



II. THREE-BODY SYSTEM WITH ZERO-RANGE 
POTENTIALS 

A. Hyper-spheric coordinates 

Let us first introduce the hyper-spheric coordinates 
{p, fii} suitable for a description of a three-body system. 
If Wj and Vi refer to the i-th particle then the hyper- 
radius p and the hyper-angle on are defined in terms of 
the Jacobi coordinates jq and y,; as [Il8| 



T = T n 



Imp- 



■A 2 , 
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T = — '"~ B/a ^ /2 
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(9) 



A 2 



sm(2aj) aa- sin 



Hi 



(di) cos 2 (ai) ' 



where l Xi and l yi are the angular momentum operators 
related to Xj and y,;. 



B. Hyper-spheric adiabatic expansion 

Let us now consider a three-body system with the 
Hamiltonian (|^). We start with the general hyper-spheric 
adiabatic expansion [^9) of the three-body wavefunction 



*(p,fi) = ^^/ n (p)$ n (p,n) 



(10) 



in terms of the complete basis $„ of the solutions of the 
hyper-angular eigenvalue equation 

\A 2 -x n ( P ) + ^J2j^ Q n *»(p,n) = o, (11) 

where Vi is the potential between particles j and k (i,j,k 
is the cyclic permutation of 1,2,3). If we truncate the 
infinite sum in ( |lo| ) we shall, according to the variational 
principle, obtain an upper bound on the discrete spec- 
trum of the system. If then the truncated expansion 
provides a collapse the full wavefunction will collapse as 
well. It is then sufficient to consider only the lowest term 
in the expansion - the so called hyper-spheric adiabatic 
approximation, 



(12) 



where the angular coordinates correspond to the " fast" 
subsystem while the hyper-radius p represents the " slow" 
subsystem. Within the adiabatic approximation the low- 
est eigenvalue A(p) of the "fast" subsystem simply serves 
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as the effective potential for the " slow" hyper-radial sub- 
system, 



d 2 A Q) + 15/4 2mE 

W + J 2 rF 



f(p) = , (13) 



where we have neglected one term, originating from the 
derivatives of the angular functions Jl^] , which is unim- 
portant for the following discussion of the collapse. 

For zero-range potentials the Faddeev equations pro- 
vide a more convenient basis for analytic insights into the 
properties of the three-body system [pjfl . The Faddeev 
decomposition of the angular wavefunction $ is 



sin(2a)$(a) = (f>(a) + —=a<t>{-) + 0(a 2 ) . (19) 
v 3 3 



Substituting ( p~9|) and (0) into (H) leads to the eigen- 
value equation |19(| for v 



-i/cosOf ) + ^sin(i/f ) _ p 1 



sin(i/^) 



(20) 



The solution v(p) of this equation defines the needed adi- 
abatic potential {y 2 — 1 /4) / p 2 for the hyper-radial equa- 
tion (|l^) from which one obtains the total energy and 
the radial wave function of the system. 



(14) 



where the three Faddeev components 4>i satisfy the three 
coupled Faddeev equations 



(A 2 - A&O) 



<f>i 



sin(2«i) 



2mp 2 , . 

_^L^$ = , i = l,2,3. (15) 



The system of Faddeev equations^ 15|) is equivalent to 
the original Schrodinger equation (Jll ) . 

C. Adiabatic solutions for the zero-range potentials 

The zero-range potentials vanish identically except at 
the origin and we are therefore left with the free Faddeev 
equations 



d 2 

da 2 ~ V ^ ' MP> a i) = ' 



(16) 



where v 2 = A + 4, and where we have restricted each of 
the components of the wavefunction to s- waves only. The 
solutions, which obey the boundary condition <fri(p, §) = 
0, are 



cj>i(p,on) = A4 sin 



(17) 



For three identical particles the index i can be dropped 
completely. For the fermionic case we imply instead that 
we have a system of, say, two neutrons and a proton in a 
state with spatially symmetric wavefunction. 

The zero-range potential appears as a boundary con- 
dition at a = [pj 



1 da<I>\ 
a$ da J 



(18) 



where p = (l/m)(mim2/(mi +ma)) is the reduced mass 
of the two particles in units of the mass m, and a is 
the scattering length. The total angular wavefunction $ 
can be expressed for small angles within a given Jacobi 
system as p9[ 



III. CORRELATIONAL COLLAPSE OF 
MANY-BODY SYSTEM 



In the small distance region, p <C a, the eigenvalue 
equation ( po|) has an imaginary root VQ = ib, where b — 
1.006, leading to an effective potential in the hyper-radial 
equation which in this region behaves as (v 2 — 1 /4) / p 2 = 



-C/p 2 , where C -- 
then becomes 



1.262. The hyper-radial equation 



dp 2 



2mE 
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f(p) = , p < a. 



(21) 



The (negative) energy E = —h 2 n 2 /(2m) is negligible 
compared to the potential when the distance is suffi- 
ciently small, p <C The hyper-radial equation then 
turns into 



dp 2 



(22) 



which has solutions of the form f(p) ~ p n , where n = 
1 ± - AC = 5 ± v . For C > j the exponent n 



, - , x , J = 5 ± v Q . For C ; 
acquires an imaginary part ±ib, i.e. 



f(p) ~ Vpe 



±ib\n(p) 



(23) 



This wavefunction has independent of energy infinitely 
many nodes at small distances or, correspondingly, in- 
finitely many lower lying states at smaller distances 
(Thomas effect). In other words, there is actually no 
finite ground state - the three-body system collapses. 

For the many-body system the energy minimum, minus 
infinity, is therefore reached at the configuration, where 
the many-body system is clusterized into separated three- 
body subsystems. The density of the many-body sys- 
tem in this configuration is then represented by a sum 
of delta-functions corresponding to collapsed three-body 
clusters. The Hartree-Fock ansatz for the wavefunction, 
however, explicitly excludes such clusterized configura- 
tions. It is therefore safe to use the Hamiltonian (Q) with 
additional potentials (H) or (R) as long as the model space 
is restricted to the Hartree-Fock product wavefunctions. 
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Let us now consider the extended Hamiltonian which 
includes the additional potential W3 either as the three- 
body zero-range potential (H) or as the density depen- 
dent zero-range potential un). Applied to a separated 
three-body cluster both potentials are non- vanishing only 
when all three constituents are located at the same point 
in space. This configuration corresponds to p = or 
ln(p) = —00. These additional potentials, therefore, will 
only change the boundary condition at p = 0. However, 
independent of the boundary condition at ln(p) = —00, 
the infinitely many nodes of the wave function (^3|) re- 
main and, therefore, the system still collapses. 

The repulsive zero-range potentials (||) or (0) need 
three particles located at the same point to provide the 
stabilizing contribution to the energy of the system. Un- 
fortunately this configuration, as we have shown, is un- 
stable under collapse. If, however, this three-body col- 
lapse is somehow removed the four-body correlations (or 
higher) will not destabilize the many-body system. In- 
deed the repulsive potentials @ or (Q) provide a contri- 
bution to the energy, which has a higher density depen- 
dence than that from the attractive two-body potential 
and hence high density configurations will not be energet- 
ically favored. The removal of the three-body collapse, 
which allows the repulsive potentials (0) or ([?]) to con- 
tribute, will therefore also stabilize the four and higher 
order correlations. 



IV. DISCUSSION 

We have shown that zero-range potentials, applied to 
a many-body system in coordinate space, lead to a spe- 
cific collapse of the system driven by three-body corre- 
lations. Neither the density dependence of the potential 
nor the three-body zero-range potential can remove this 
correlational collapse. The Hamiltonian (En, perhaps ex- 
tended by the additional potential ^ or (Q), can still 
produce meaningful results, provided the allowed vari- 
ational model space is chosen consistently. Any three- 
body correlation must be excluded a priori from the 
wavefunction of the many-body system approximated by 
the Hamiltonian (jsj) . Going beyond the Hartree-Fock 
(anti)symmetrized product wavefunction would immedi- 
ately introduce dangerous effects of three-body correla- 
tions. 

To remove the three-body collapse one has to regular- 
ize the zero-range potential by introducing a finite length 
scale R, which will alter the —Cj p 2 behaviour of the effec- 
tive potential at p ~ R. Such regularization is automat- 
ically achieved by finite-range two-body potentials. The 
three-body system will then have approximately ln(a/i?) 
bound states with the ground state having a finite bind- 
ing energy of the order of ft 2 / (2mR 2 ). 

When the three-body system is regularized the corre- 
lations of higher number of particles can be safely in- 



cluded in the many-body wavefunction as the ordinary 
additional potentials (0) or (^) will make the system sta- 
ble. 
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